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ABSTRACT 

We present dynamically consistent solutions for hot accretion onto unmagnetized, rotating white dwarfs 
(WDs) in five quiescent dwarf novae. The measured WD rotation rates (and other system parameters) 
in RX And, SS Cyg, U Gem, VW Hyi and WZ Sge imply spindown of the WD by an extended hot flow 
emitting most of its X-rays in the vicinity of the stellar surface. In general, energy advection is absent 
and the flow is stable to convection and hydrodynamical outflows. In rapidly rotating systems, the X-ray 
luminosity provides only an upper limit on the quiescent accretion rate because of substantial stellar 
spindown luminosity. We suggest that the presence of hot flows in quiescent dwarf novae may limit the 
long-term WD rotation rates to significantly sub-Keplerian values. 

Subject headings: X-ray: stars - binaries: close - accretion, accretion disks - stars: white dwarfs 


1. INTRODUCTION 

Dwarf novae (DN) are accreting binary star systems 
with a white dwarf (WD) primary fed by a main-sequence 
donor via Roche-lobe overflow. They represent a sub¬ 
class of cataclysmic variables (CVs) and experience semi¬ 
regular, luminous outbursts, during which accretion onto 
the WD proceeds at a high rate. Most of the time, how¬ 
ever, DN are in quiescence — a phase during which the 
accretion rate onto the WD (and the system luminosity) 
is considerably reduced (see Warner 1995 for a review). 
According to the disk instability model (DIM), such a be¬ 
havior arises because most of the mass transferred by the 
companion star builds up in an unsteady disk during qui¬ 
escence (Cannizzo 1993; Lasota 2001). 

Quiescent DN are hard X-ray sources with typical lu¬ 
minosities of ~ 10^° — 10^^ erg s“^ (see, e.g., Cordova & 
Mason 1983; Patterson & Raymond 1985). Spectral fits to 
this X-ray emission suggest a Bremsstrahlung origin from 
gas with temperatures ~ 2 — 20 keV (Patterson & Ray¬ 
mond 1985; Eracleous et al. 1991; Belloni et al. 1991; 
Yoshida et al. 1992; Mukai & Shiokawa 1993). This X- 
ray emission is commonly attributed to the boundary layer 
(BL) at the interface between the WD and the thin accre¬ 
tion disk around it. At low accretion rates ( ^ 10^® g s“^, 
typical of quiescent DN), the gas in the BL is hot and op¬ 
tically thin, hence it is a substantial source of hard X-ray 
emission (Pringle & Savonije 1979; Tylenda 1981; Patter¬ 
son & Raymond 1985). Detailed calculations by Narayan 
& Popham (1993) show that the optically-thin BLs of disk- 
accreting WDs can also be radially extended (of order a 
WD radius) and that energy advection is an important 
element of their internal structure. 

On the other hand, Meyer & Meyer-Hofmeister (1994) 
presented two arguments in favor of a more extended, in¬ 
ner hot flow structure in quiescent DN. They noted that: 
(i) the standard DIM predicts accretion rates at the WD 


surface which are smaller than those inferred from quies¬ 
cent X-ray luminosities by typically more than one order 
of magnitude and (ii) the presence of an extended, low- 
density hot flow in the WD vicinity may explain the ob¬ 
served delay (of ~ 0.5 — 1 day) in the rise to outburst of 
the EUV light relative to the optical light (e.g., Mauche, 
Mattei & Bateson 2001, but see Smak 1998). Our current 
understanding of the structure and properties of quiescent 
disks is rather limited (Menou 2002), so studying the pos¬ 
sible presence of an extended, hot flow in quiescent DN is 
certainly worthwhile. 

Recently, Medvedev & Narayan (2001a) discovered solu¬ 
tions for hot accretion onto unmagnetized, rotating com¬ 
pact stars. They found that at large stellar spin rates, 
dissipation in the hot flow is dominated by stellar brak¬ 
ing in a “hot settling flow” (HSF) configuration, whereas 
in the opposite limit, the flow reduces to a conventional 
advection-dominated accretion flow (ADAF; Narayan & 
Yi 1994; 1995). Both flows exist at relatively low accre¬ 
tion rates (^ 10“^ of Eddington). Medvedev & Narayan 
(2001b) have also addressed the thermal stability of their 
(cooling-dominated) settling solution, and found that it 
is most likely stable, once turbulent thermal conduction 
is accounted for. The work of Medvedev & Narayan was 
largely focused on accretion onto neutron stars, but the 
HSF solutions are also valid for WD accretion. This is 
important because, in general, much more is known about 
WDs in CVs than about neutron stars in close binary sys¬ 
tems. In particular, the rotation rates of several WDs in 
DN have been measured during the last few years thanks 
to HST spectroscopy (see Sion 1999 for a review). In this 
Letter, we apply the model of Medvedev & Narayan to 
WDs in quiescent DN by constructing numerical solutions 
for hot accretion in five systems with relatively well known 
system parameters. 

2. SYSTEM PARAMETERS 
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We have gathered existing data on quiescent DN for 
which estimates of the WD rotation rate and the quies¬ 
cent X-ray luminosity were both available. A summary 
of the system parameters adopted for this study is given 
in Table |^. For the WD masses, M^d, and orbital in¬ 
clinations, z, we uniformly adopt the values given in the 
catalog of Ritter & Kolb (1998). To derive the WD ro¬ 
tation rates, we use the observationally inferred values of 
Urotsini compiled by Sion (1999) and make the assump¬ 
tion that the stellar rotation axis is aligned with the ro¬ 
tation axis of the orbital motion. A spin parameter, s, 
can then be calculated as the ratio of Uot to the Keple- 
rian rotation rate at the WD radius, VK(Rwd)- For the 
WD radii, we simply use the reasonable mass-radius re¬ 
lation: i?wd = 5 X 10® cm (Mwd/1.2 Mq)”^ ®®. Finally, 
references from which quiescent X-ray luminosities were 
collected can be found in Table |^. When several values 
existed for a given system, we generally adopted the value 
corresponding to the hardest X-ray spectral range (given 
the typically hard Bremsstrahlung spectrum of quiescent 
DN). Note that the accuracy of the system parameters 
adopted is not crucial because we only wish to illustrate 
the general properties of hot accretion in these systems.^ 

3. HOT ACCRETION SOLUTIONS 

To investigate the nature of hot accretion in quies¬ 
cent DN, we use the numerical code described in detail 
by Medvedev & Narayan (2001a). We solve the height- 
integrated, axisymmetric hydrodynamical equations for 
the conservation of mass, radial and angular momenta, and 
two energy equations for the proton and electron fluids. 
These equations are solved by the relaxation method with 
appropriate boundary conditions. A special grid was used 
to carefully resolve steep gradients near the stellar surface. 
At the outer boundary, the electron and proton tempera¬ 
tures and the angular velocity of the flow are fixed to the 
values appropriate for an ADAF. At the inner boundary, 
where the flow meets the stellar surface, the flow tempera¬ 
ture and the angular velocity are forced to match those of 
the star. We assume a stellar temperature of ^ few x 10® K 
(the exact value is unimportant here). The mass accretion 
rate, M, is a free parameter, and a standard Shakura- 
Sunyaev a-viscosity is used. We assume that most of the 
viscously generated heat goes to protons. However, this 
is unimportant because the electrons are coupled to the 
protons via Coulomb collisions very efficiently. The elec¬ 
trons are cooled via optically thin Bremsstrahlung emis¬ 
sion. Thermal conduction and (self-absorbed) synchrotron 
cooling are not included, for simplicity (the latter should 
be negligible). The effects of Comptonization of the stellar 
radiation and self-Comptonization of the Bremsstrahlung 
radiation are included, but both are very weak. We gen¬ 
erally adopt an adiabatic index 7 = 1 . 6 , close to the 
7 = 5/3 of a non-relativistic, monoatomic ideal gas. In 
all models we choose the outer radius of the flow to be at 
~ few X 10®® cm to better isolate the effects of the inner 
and outer boundary conditions. In real systems, the outer 
radius is probably much smaller. The adopted parame¬ 
ters, namely the stellar mass, M^d, its radius, i?wd, and 
its spin parameter, s, are given in Table In all systems 


but WZ Sge, we used a conventional value of the viscos¬ 
ity parameter, a = 0.1. In WZ Sge we set a = 0.02, as 
discussed below. 

Figure ^ shows the numerical solutions obtained for the 
five systems listed in Table |^. In each case, the accretion 
rate in the steady, hot flow has been adjusted so that the 
total Bremsstrahlung emission fits the quiescent X-ray lu¬ 
minosity given in Table |^. The four panels in Figure 
show the radial profiles of density, p (in g cm“®), elec¬ 
tron and proton temperatures, Te and Tp respectively (in 
Kelvin; the curves overlap because « Tp indicating ef¬ 
ficient Coulomb transfer of energy between the ions and 
the electrons), angular rotation velocity, Q (in units of the 
Keplerian value at the WD radius, flK{Rwd)) and radial 
velocity, Vr (in units of the speed of light, c). The profiles 
are approximate power laws, with narrow boundary layers 
close to the WD surface (except for D, which smoothly 
matches the stellar rotation rate). The profiles of fl and 
Vr indicate that boundary effects are present over a sub¬ 
stantial range of radii, at both the inner and outer edges. 

The scalings in an ADAF are p oc i?“®/^ and Vr oc 
/j-i /2 (Narayan & Yi 1994), while in a HSF they are 
p oc R~^ and Vr oc i?® (Medvedev & Narayan 2001a). 
The profiles of p and Vr in Fig. ^ deviate from perfect 
power laws and are intermediate between the above two 
solutions. The Bernoulli number is negative everywhere 
(for 7 ~ 5/3), except in a narrow region at the outer edge 
of the flow where the ADAF conditions hold (in an ADAF, 
Be > 0). Hence, the solutions should be convectively sta¬ 
ble (Narayan et al. 2000; Medvedev & Narayan 2001a). 

The solutions shown in Fig. ^ are not “advection- 
dominated” in the sense that the gas entropy decreases 
inward (it increases inward in an ADAF). The flow en¬ 
ergetics is best described by computing the ratio of the 
integrated Bremsstrahlung luminosity to the available ac¬ 
cretion luminosity, Lacc = GM^^M/R^^d- The accretion 
rates for the five systems are quite low (Table ||). The 
corresponding ratios, rjes = Tx/Tacc, are ~ 1 for low-s 
systems and larger than unity for large s, indicating sig¬ 
nificant extraction of rotational energy from the central 
star. Thus, in large-s systems such as WZ Sge, the M 
value given in Table should be considered as a rough 
upper limit because Lx is essentially powered by stellar 
spindown with little contribution from the accretion lumi¬ 
nosity. However, the luminosity of a HSF is sensitive to a. 
We were able to fit the observed Lx of WZ Sge only for 
a « 0.02. Perhaps, this may indicate that the value of a 
may be lower than generally expected. 

The value of the adiabatic index has a strong effect 
on the flow structure. We illustrate this by calculating 
a model for U Gem with 7 = 4/3 (all other parameters 
being the same). It is shown in Fig. ^as a thick solid line. 
At large radii, the flow structure now adopts an ADAF 
configuration, with p oc i?“®/^, Vr oc R~^G and a larger, 
constant value of LljLlK- The Bernoulli number is positive 
everywhere in the hot flow (excluding the BL) and en¬ 
ergy advection is dominant in the outer regions of the flow 
with ADAF properties. However all the advected energy 
is now radiated in the optically thin BL so that Tyeg « 1 
as well. This influence of 7 is consistent with the results 


^For instance, even if the WD in WZ Sge was much more massive than assumed in Table |l| (e.g. Steeghs et al. 2001), it would still be the 
fastest rotator among all five systems, as measured by the s parameter. 
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of Medvedev & Narayan (2001a). According to their clas¬ 
sification, the solutions we obtained for 7 = 1.6 and the 
inner zone of the flow in the case with 7 = 4/3 are “settling 
ADAFs” (WZ Sge clearly corresponds to a HSF, however). 

Fig. ^ shows the Bremsstrahlung luminosity, per unit 
logarithmic fractional distance from the WD surface, for 
the five solutions with 7 = 1.6 shown in Fig. |^. The den¬ 
sity in the flow is low enough that Comptonization has a 
negligible effect on the emitted spectrum, even in the dens¬ 
est regions, near the WD surface. Fig. || clearly shows that 
X-ray emission originates mostly from regions in the close 
vicinity of the WD (the same is true for the U Gem model 
with 7 = 4/3). The brightest regions of the flow are also 
approximately the hottest, with temperatures ~ 10® K for 
WZ Sge, ~ 2 X 10® K for VW Hyi and ~ 4 x 10® K for the 
three other systems. 

4. DISCUSSION 

We constructed numerical models for hot accretion onto 
unmagnetized, rotating WDs in five quiescent DN. Our so¬ 
lutions are a significant improvement over previous work 
on hot accretion in this context (Katz 1977; Kylafis & 
Lamb 1982; Mahasena & Osaki 1999; Menou 2000) in that 
both the WD rotation and the viscous nature of the flow 
are accounted for in the present case. This makes these 
solutions plausible modes of accretion in quiescent DN. 

The outer boundary conditions chosen are somewhat 
artificial and constitute a significant shortcoming of the 
present study. It is unclear exactly how the transition 
from an inner hot flow to an outer thin disk (known to 
be present) would occur. Strong boundary effects may 
be expected for a hot flow with a limited radial extent. 
Nonetheless, our work has the virtue of isolating the main 
properties of the hot flow, independently of the disk. It 
will be important in the future to study the structure of 
a hot flow with a smaller radial extent and more realistic 
boundary conditions. For the non-relativistic gas consid¬ 
ered here and non-dominant magnetic fields, we expect 
rather large values of 7 and therefore flows which are not 
subject to outflows (Blandford & Begelman 1999) or con¬ 
vection (Narayan et al. 2000; Quataert & Gruzinov 2000). 

Menou (2000) proposed that energy advection in a hot 
flow could be important in powering a dominant EUV 
emission component that could explain the strong He II 
emission lines observed in many quiescent DN. Since the 
hot flow solutions presented here lack energy advection 
(for 7 ’s near 5/3), this possibility seems unlikely. 

While a star is spun up by accretion via a disk bound¬ 
ary layer (except when rotating near breakup; Popham & 
Narayan 1991; Paczynski 1991), it is spun down when ac¬ 
creting via a hot flow like an ADAF, a HSF or the hot flow 


solutions presented here (except at very low spin rates; 
Medvedev & Narayan 2001a). This property has poten¬ 
tially important consequences for the spin history of WDs 
in DN. Sion (1999) emphasized that in the standard pic¬ 
ture WDs should rotate near break up in DN, but observa¬ 
tions suggest that rapidly spinning WDs are rare. If a WD 
is spun down via hot accretion during quiescence, an equi¬ 
librium at substantially sub-Keplerian rotation rates may 
be expected from the balance between spin-up in outburst 
and spin-down in quiescence. This possibility requires fur¬ 
ther study with more appropriate boundary conditions. 
We note that other mechanisms could contribute to spin¬ 
ning down single WDs (Spruit 1998) and WDs in DN (and 
thus get rid of the angular momentum acquired during ac¬ 
cretion), including angular momentum losses induced by 
coupling to an expanding envelope during nova explosions 
(Livio & Pringle 1998; Sion et al. 2001). 

X-ray eclipses in several quiescent DN indicate that the 
emission originates from the vicinity of the WD (Wood et 
al. 1995; Mukai et al. 1997; Pratt et al. 1999; Ramsay 
et al. 2001). Without a detailed comparison, it is unclear 
whether the hot flow solutions presented here satisfy the 
existing observational constraints or not. Such a compar¬ 
ison should take into account the spectral coverage of the 
various X-ray instruments used and it probably requires 
a reliable boundary layer model taking into account the 
role of heat conduction in potentially modifying its emis¬ 
sion properties (which is beyond the scope of the present 
paper). It is certainly encouraging that X-ray emission is 
highly concentrated near the stellar surface in our models, 
but it is presently unclear whether these models will be 
able to account for the available X-ray eclipse data or not. 
We note, for instance, that the X-ray eclipse of OY Gar 
observed by Ramsay et al. (2001) suggest an X-ray flux 
approaching zero at eclipse minimum and a small verti¬ 
cal extent for the X-ray emitting region. These properties 
may not be easily reconciled with the hot flow models pre¬ 
sented here. In the future, it will also be possible to probe 
the structure of hot flows in (even non-eclipsing) quies¬ 
cent DN with detailed X-ray spectroscopic diagnostics, as 
illustrated in Menou, Perna & Raymond (2001; see also 
Narayan & Raymond 1999). 
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Table 1 

SYSTEM PARAMETERS 


System 

name 

(Me) 

/?wd 

(10® cm) 

i 

n 

Cot 

(km s“^) 

5 

(Cot/kk,wd) 

Lx 

(erg s“^) 

M 

(g s-i) 

^eff 
i -^acc) 

RX And 

1.14 

5.1 

51° 

190 

0.03 

3 X lO^i “ 


1.0 

SS Cyg 

1.19 

5.0 

37° 

500 

0.09 

1032 b 

3.3 X 10^4 

1.0 

U Gem 

1.26 

4.9 

70° 

< 110 

< 0.02 

1031 c 

3 X 1013 

1.0 

VW Hyi 

0.63 

6.2 

60° 

460 

0.13 

7 X lO^o ^ 

3.5 X 1013 

1.5 

WZ Sge 

0.45 

7.0 

75° 

1240 

0.42 

3 X 10^0 « 

1.1 X 1013 

3.2 


NOTE. - (a) Eracleous, Halpern & Patterson (1991) (b) Yoshida, Inoue & Osaki (1992); Mukai & Shiokawa (1993) (c) 
Szkody et al. (1996) (d) Mukai & Shiokawa (1993); Eracleous et al. (1991) (e) Eracleous et al. (1991); Mukai & Shiokawa 
(1993). 
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log D = log[(R-R^J/R^n] 


Fig. 2.— Bremsstrahlung lumimpsity of the hot flow, per unit logarithmic fractional distance from the white dwarf surface, for the five 
models with 7 = 1.6 shown in Fig m (same notation). Most of the X-ray emission originates from the white dwarf vicinity. 





